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Relativistic heavy-ion collisions lead to a final state which has a higher degree of strangeness 
saturation than those of elementary collisions. A systematic analysis of this phenomenon, based 
on the strangeness saturation factor, 7^, is made for C+C, Si+Si and Pb+Pb collisions at the 
CNj ' CERN SPS collider and for Au+Au collisions at RHIC energies. Strangeness saturation is shown to 

I increase with the number of participants within a colliding system, at both CERN SPS and RHIC 

. energies. The saturation observed in central collisions of lighter nuclei deviates from that seen in 

(N ■ peripheral collisions of heavier nuclei with an equivalent participant number, which could be due to 

^ ' the difference in nuclear density. 
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I. INTRODUCTION 



It has been shown that statistical-thermal models are able to reproduce the multiplicities measured in relativistic 
heavy-ion collisions with remarkable success. This is accomplished with a very small number of parameters- the 
. temperature, baryon-chemical potential, fiB, and a factor measuring the degree of strangeness saturation, 7^. As is 
' now well known, there is very little difference between the temperatures observed in p + p and relativistic heavy-ion 
collisions. The extracted strangeness saturation factor, 7^, is however very different in p + p and heavy-ion collisions. 
In this paper we focus on hadron multiplicities and extract the thermal parameters as a function of system-size and 
T— I . energy. 

C\| • The recent study in |^ (cf. table III therein) irnpressively demonstrated that, with increasing system-size at SPS 
[ energies, the strangeness saturation increases. In [Q| we have shown that at a beam energy of 158 AGeV, in collisions 
t of lead-on-lead nuclei, the strangeness saturation continuously increases with centrality. However, strangeness (as 
rS I measured by fully-integrated kaon and antikaon multiplicities) is clearly below saturation. A preliminary analysis 
' 1^7 @i D of the centrality dependence at RHIC energy of \/^jvAf ~ ^"^^ GeV points to a further increase of strangeness 

towards saturation for central collisions of gold nuclei. An independent analysis Q confirms this finding. 
(— I \ This paper is divided into several sections. Firstly the system-size dependence of the thermal parameters is de- 
■ termined using 47r-yields from central C-l-C and Si-|-Si collisions and centrality-binned Pb-|-Pb collisions [|[ |l^ 
.'^ l at 158 AGeV at the CERN SPS. For comparison, centrality-binned mid-rapidity yields from Au+Au collisions at 
^ ' V^NN — GeV and Pb+Pb coUisions at SPS energy |^ are analysed, despite the danger in applying the 
5-H thermal model to yields in a limited rapidity window. Finally, the energy dependence of the thermal parameters is 
. 5^ , further elucidated by analysis of central Pb+Pb yields measured by NA49 at 40, 80 and 158 AGeV |l|, |l3[ |l|, |l|, ^ . 



II. ANALYSES OF HADRON MULTIPLICITIES 

In the thermal model, hadron multiplicities can be described ||l^, |l^, |l^ |2^, ^ by the grand-canonical partition 

function Z{V,T,fli) = Tr{e t }, where H is the statistical operator of the system, T denotes the temperature, 
and Hi and Qi represent the chemical potentials and corresponding conserved charges respectively. In the analysis 
of 47r-data, the net-zero strangeness and the baryon-to-electric charge ratio of the colliding nuclei constrain the 
components of /ii = {^b, fJ'S, ^J■Q)■ The particle numbers are given, in the Boltzmann approximation, by 

iVr" = ^5.71^-1 / dm, e-^^BW(m,), (2.1) 
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where we include phenomenologically a strangeness saturation factor, 7^, with \Si\ the number of valence strange 
quarks and anti-quarks in species i |22 (i.e. 7s for the kaons and 7^ for cj)) to account for incomplete equilibration in 
this sector, Ei = \/p^ + rnf, and BW is the Breit-Wigner distribution. The particle numbers to be compared with 
experiment are Ni — Nf""^ + Br(j/ i)N^""^, due to decays of unstable particles with branching ratios Br(j 



For small participant numbers (typically A'part below 40), one has to resort to a canonical or micro-canonical formalism 

MM. 



A. System-size dependence 

1. Analysis of fully-integrated yields 

In order to extract the system-size dependence of the thermal parameters we analyse 47r-multiplicities of tt^ , , 
(f> and iVpart (taken as the sum over all baryons) in 6 centrality bins in the reaction Pb-|-Pb Jp| Jl(|_ (at 158 AGeV) 
and for central Si-|-Si and C-l-C collisions [|| at the same energy. Our previous analyses [|[ |^|j7|ll of the Pb-|-Pb 
system included p yields. They are excluded in this analysis in order that the Pb-, C- and Si systems be treated 
equivalently. No weak feed-down corrections have yet been applied to the peripheral Pb-, or C- and Si systems [p5| . 
Due to the rather limited data set, the freeze-out temperature was fixed at 165 MeV, independent of centrality and 
colliding system. This is supported by a variety of fits to both heavy-ion and elementary collision systems in this 
energy regime ^ ^ Owing to the size of the C- , Si- and peripheral Pb systems, strangeness was treated 

canonically in all systems. As shown in for systems of this size at SPS energy, it is sufficient to treat the baryon- 
and charge content grand-canonically. The results are displayed in Figs. ^ and |^ with the specifics of each fit explained 
in the captions. 

The strangeness saturation factor, 7^, shows an increasing trend with collision centrality in the Pb-|-Pb system, 
except possibly over the two most central bins (see Fig. |l]). It is also clear that the C-l-C and Si-|-Si systems lie above 
the trend suggested by the Pb-|-Pb points. This suggests that peripheral Pb-|-Pb coUisions are not equivalent, with 
respect to strangeness saturation, to central collisions of lighter nuclei with the same participant number. In the C-f-C 
and Si-|-Si systems the baryon chemical potential is also lower than in the peripheral Pb-|-Pb bins (refer to Fig. ||). It 
should be stressed that the only direct baryon information we have at our disposal in this analysis is the number of 
participants. It would appear that, in the Pb-f Pb system, /is decreases as the collisions become more central. As can 
be seen in Fig. the inclusion of p yields in the Pb analysis (squares) leads to a roughly centrality-independent baryon 
chemical potential of approximately 250 MeV. In order to investigate its interplay with 7^, we fixed /is at 250 MeV 
in all collisions. This affected appreciably only the most central Pb-|-Pb bin, resulting in a monotonic increase in 7s 
within the Pb system. Thus, the drop in 7s over the most central bins of the Pb system is driven by iVpart which, 
due to its small error relative to other species in the most central bin, is weighted heavily in a x^-analysis. Since the 
errors, particularly in the C and Si systems have not yet been well-established we repeated the fits minimising 
the 'quadratic deviation' [|9| defined by: 

V {Mr''") 

where M^^^ and M™°'*^' are the experimental and model-predicted multiplicities of hadron species i respectively. 
The results of this analysis are shown in Figs. |l| and ^ as triangles. Again 7s in the Pb+Pb system shows a monotonic 
trend, while the C+C and Si+Si systems still show a sizeable deviation from the peripheral Pb points. This 
cannot be attributed to rescattering, since it is more dominant in peripheral Pb-|-Pb than in central C-l-C and Si+Si 
reactions p^ . In |8| it is shown that the strangeness enhancement, as measured by the ratios of strange to non-strange 
mesons, in these systems scales with /2, the fraction of participants which undergo multiple collisions. As shown in 
Fig. H, 7s scales with this variable too. The strangeness saturation extracted from p-l-p collisions at ^/s = 19.4 GeV |2^ 
(denoted in Fig. ^ by the square) suggests a strong flattening off of 7s for small systems. What is surprising is 
the approximate equivalence of /2 and 7s for a number of points in Fig. |l| (/2 in this figure is denoted by the diamonds). 



2. Mid-rapidity analysis 

When applying the thermal model to 47r-data, many dynamical effects cancel out in ratios of the fully-integrated 
hadron yields [E9|. In particular, effects due to ffow disappear if the freeze-out surface is characterized by a single 
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temperature and chemical potential. When applying the thermal model to mid-rapidity data this is true only when the 
Bjorken model ]30[| holds. Furthermore, in a limited rapidity window there is no guarantee that the total strangeness 
should be zero, nor that the baryon-to-charge ratio in this kinematic region should be set by the colliding nuclei. In 
a grand canonical approach this affects the constraints on the chemical potentials. In the canonical formulation the 
'canonical suppression volume' (i.e. the volume in which quantum numbers are exactly conserved and that used to 
calculate the densities) and the 'normalisation volume' (i.e. the volume required to convert densities to yields) are 
not necessarily the same. A further complication arises in the treatment of decays. The experimental yields include 
feed-down from heavier resonances into stable, final-state particles. With mid-rapidity data, this requires careful 
consideration of the decay kinematics, as particles in a certain rapidity range will in general decay into particles in 
different kinematic windows. 

Despite these disclaimers we analyse the following mid-rapidity yields measured by the NA49 and PHENIX collab- 
orations: 

i) NA49 mid-rapidity yields of tt^, if^, p and p in the reaction Pb(158 AGeV) + Pb in 6 centrality bins 

ii) PHENIX mid-rapidity densities of tt*, and p^ in the reaction Au + Au at = 130 AGeV in 5 centrality 
bins ||ll|. 

The PHENIX yields were not corrected for weak decays. PHENIX estimate the probability for reconstructing 
protons from A decays as prompt protons at 32% aX, pT — ^ GeV/c [|ll|. The PHENIX analysis was performed with 
both 0% and 50% feed-down from weak decays, while no weak feeding was included in the SPS analysis. In all cases 
the grand-canonical formalism was applied, with the total strangeness set to zero but pq fit as a free parameter. 

In Fig. 1^ the system-size dependence of 7s at mid-rapidity is shown for SPS (left panel) and RHIC (right panel). 
In the SPS plot the results of our earlier analysis of fully-integrated NA49 yields [|| are included for comparison. 
It is seen that 7^, as extracted from the mid-rapidity NA49 data, is well above that obtained from the analysis of 
the fully-integrated NA49 yields. In order to exclude the possibility that the difference in the strangeness saturation 
extracted from the 47r- and mid-rapidity data analysed here is due solely to different strange hadrons included in the 
fits, the 47r NA49 analysis was repeated with the hidden-strangeness (j) excluded. In this way the two NA49 analyses 
are equivalent with respect to strange particles. This led to a slight decrease in 7^ in the most peripheral bins, and 
thus an even larger difference between mid-rapidity and fully-integrated results. Thus, certainly at SPS energies, the 
degree of strangeness saturation is far higher in the central rapidity region. Included in the RHIC plot of Fig. | is /a. 
This fraction of multiply-struck participants parametrises the system-size dependence of the strangeness saturation 
factor in Au-|-Au collisions at RHIC energy remarkably well. 



B. Energy dependence 

In order to further investigate the energy dependence of the thermal parameters we analyse the fully-integrated 
yields of tt^, K^, A and A in central Pb+Pb collisions at 40, 80 and 158 AGeV, supplemented with K^, S~, "E^ and 
(/) multiplicities at 158 AGeV ^ |l|, |l|, 0. In Figs. | and |, 7^ and the Wroblewski factor [|l|' which measures 
the ratio of newly created ss pairs to newly created non-strange valence quark pairs at the primary hadron level: 

{uu) + {ad) 

are displayed as a function of the collision energy. Included in the figures are the results of our earlier system-size 
analysis of the Pb-|-Pb system at CERN SPS and the Au-|-Au system at RHIC [1). It should be noted that the 
centrality cuts on the most central Pb-fPb collisions are slightly different at the various energies (7.2% at 40 and 
80 AGeV, and 5% at 158 AGeV). In view of the system-size dependence extracted in the previous section, this will 
raise the 40 and 80 AGeV points relative to the 158 AGeV point. In Fig. ^ one observes that As for central collisions 
decreases with collision energy from 40 AGeV (in agreement with ^^). Within a given collision system, it shows a 
systematic increase with participant number, while remaining above the typical value of 0.2 seen in pp collisions [^ . 
With respect to 7s there is fairly good agreement between the most peripheral heavy-ion bins and the results from 
elementary systems of comparable energy . 

For comparison we also show the value of 7s extracted from mid-rapidity Pb+Pb data ||l^, ^ |l^ at 158 AGeV 
(open circle) in Fig. ^. As can be seen, the strangeness saturation in the mid-rapidity region is clearly greater than 
that averaged over 47r. 
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III. SUMMARY 



In conclusion, the strangeness saturation factor, 7s, has been shown to increase with participant number in the 
Pb+Pb system at the CERN SPS as well as the Au+Au system at RHIC. Central coUisions of C+C and Si+Si at 
SPS energies deviate, with respect to strangeness saturation, from peripheral Pb+Pb collisions. However, 7^ is seen 
to scale with the fraction of multiply-struck participants, /2. In fact, /2 remarkably tracks the A'part-dependence of 
7s as extracted from mid-rapidity yields in Au+Au collisions at RHIC. Where both mid-rapidity and fully- integrated 
data was available, the degree of strangeness saturation observed at mid-rapidity was found to be consistently higher 
than that extracted from 47r-data. 
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FIG. 1: The system-size dependence of the strangeness saturation factor, 7s, as extracted from centrality-binned Pb+Pb |^ 
and central C+C and Si+Si data under various fit conditions. The circles with error bars represent the results of our x^- 
analysis assuming 50% feeding from weak decays, while the triangles show the results minimising the quadratic deviation 
(again assuming 50% feeding). For comparison, the results of our earlier Pb system analysis jsj, with p's included in the fit, 
are included (squares). Also shown are the fraction of participants which underwent multiple collisions, /2, as extracted from 
a Glauber calculation (diamonds). 
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FIG. 2: The system-size dependence of the baryon chemical potential, /is, as extracted from centrality-binned Pb+Pb j^, p^ , 
and central C+C and Si+Si data Q under various fit conditions. The circles with error bars represent the results of our 
X^-analysis assuming 50% feeding from weak decays, while the triangles show the results minimising the quadratic deviation 
(again assuming 50% feeding) 
are included (squares) . 



For comparison, the results of our earlier Pb system analysis jEj, with p's included in the fit, 
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FIG. 3: The strangeness saturation factor, 7s, as extracted from centrality-binned Pb+Pb j^, |l^ (triangles) and central C+C 
and Si+Si data fel (circles), as a function of /2, the fraction of multiply-struck participants. The results shown are those 
obtained with /is fixed at 250 MeV and T at 165 MeV, assuming 50% weak feed-down. For comparison, the strangeness 
saturation as extracted from p-|-p collisions at ^/s = 19.4 GeV |^ is included (square). 
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FIG. 4: Left Panel: Comparison of the strangeness saturation factor, 7s, extracted from mid-rapidity NA49 data ^ (up 
triangles) with the results of our earlier analysis of NA49 47r-yields (squares) [^. Right Panel: The strangeness saturation 
observed in Au-|-Au collisions as extracted from PHENIX data jl^. The analysis was performed assuming 50% weak feed- 
down (down triangles) and 0% weak feed-down (up triangles). Also shown are the fraction of multiply-struck participants, /2, 
obtained from our Glauber calculation (dashed line). 
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158 AGeV |l^ [l^ (up triangles), together with the results of our earlier analysis of centrality-binned Pb+Pb collisions 
at 158 AGeV (down triangles) and Au+Au collisions at RHIC (diamonds) |^. For comparison we show the results obtained 
from pp colhsions (filled squares) and pp collisions (open squares) at various energies M]. The open circle is extracted from 
mid-rapidity yields [||, [|, |l|. 
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FIG. 6: The ener gy dependence of the Wroblewski factor, As, extracted from central Pb+Pb coUisions at 40, 80 and 
158 AGeV |14[, (up triangles), together with the results of our earlier analysis of centrality-binned Pb+Pb col- 

lisions at 158 AGeV(down triangles) and Au+Au collisions at RHIC (diamonds) Q]. For reference we show the typical value 
of 0.2 extracted from pp systems 0]. 



